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& ' Abstract 

We apply a recent adaptation of White's density matrix renormalisation group (DMRG) 
method to a simple quantum spin model, the dimerised XY chain, in order to assess the applica- 
bility of the DMRG to quantum systems at non-zero temperature. We find that very reasonable 
[ results can be obtained for the thermodynamic functions down to low temperatures using a very 

small basis set. Low temperature results are found to be most accurate in the case when there is 
a substantial energy gap. 

X 

Since it's recent inception, White's denstiy matrix renormalisation group (DMRG) method has 
been established as the method of choice for determining static, low energy properties of one-dimensional 
quantum lattice systems Q-Q]. Extensions to the calculation of dynamical properaties || and even to 
the study of low temperature properties of two dimensional systems [|| have been forthcoming. More- 
over, Nishino's formulation of the DMRG for two dimensional classical systems || has paved the way 
for the study of one dimensional quantum systems at non-zero temperature. In this letter we present 
what is, to the best of our knowledge, the first application of the DMRG to the thermodynamics of a 
quantum system. 

The system that we consider is a simple spin chain model — the dimerised, S — 1/2, XY model 

N/2 
i=l 

where Si is a spin-1/2 operator for site i on an (even) chain of N sites, with periodic boundary 
conditions. 







This model is exactly solvable Q ; the Helmholtz free energy tp is given by 
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where Zn = Tre ^ is the partition function, (3 = 1/T is the inverse temperature, <fi(9) = cos£(#) 
7cos(0+£(l?)), £(0) = -tan- ; 
such as the internal energy u 



7 cos(0+^(6*)), = — tan 1 ^ sl " o ^ e , 7 = J2I J\ and we have set Ji = 1. Thermodynamic properties 



^Jj$- and the specific heat cy = — (3 2 §^ are readily obtainable from 



Because of its solvability and the fact that it possesses both gapless (7 = 1) and gapped (7^1) 
excitations, (|l|) presents a useful test model for the extension of the DMRG to quantum systems at 
non-zero temperature. 

Now Nishino's formulation of the DMRG applies directly to classical, two dimensional spin systems 
and so in order to apply it to the quantum system (Q), we must first invoke the Trotter-Suzuki method 
@. That is, we make the decomposition H = H\ + H2 where 
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and hi t j = SfS/j + SfSj. We then apply the formula 
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Inserting a collection a of 2M complete sets of states into (|4|) then yields 
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where periodic boundary conditions of M+1 = a}, cr J N+1 = a{ are assumed. 

Now, because the N/2 terms in the sums (^|) commute and act on different pairs of sites, we have 
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for i — 1,2. Moreover, 
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where T = T\T2 is the virtual transfer matrix and 
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The matrices 71 and T2 are depicted graphically in Fig. 1. 
It follows that 

V> = ^ (M) = lim - l0gAmax (12) 
v v m^oo 2/3 v ; 

where A max is the eigenvalue of T with maximal modulus and generally depends on the Trotter 
dimension M. 

We may now apply the formulation of the DMRG for transfer matrices to the calculation of A max . 
We commence by defining an initial transfer matrix % for a single-site system block, which connects 
it to adjacent sites to the left and the right viz 

T a {a'n'n'\an^a"n") = J2 T i(o-'o-"\n'n")T 2 (n"n\n"fi) (13) 

n" 

The initial block, having just one site, has to = 2 states, n — 1 or 2 (f or J.) (see Fig. 2). 

We next define an initial environment block and associated transfer matrix %. in precisely the same 
way: T c = T s . We define a superblock using the system and environment blocks in addition to two 
added sites, arranged in a periodic fashion (see Fig. 3). That is, the superblock transfer matrix is 
given by 

T {(j l 1 n' l a l 2 n' 2 \ainia 2 n 2 ) = {TiT2){a' l n' 1 a 2 n' 2 \ainia 2 n 2 ) (14) 

= 7^ (crj n.' x erg I <ti m (72 ; a'[a 2 )%(a l 2 n' 2 a' 1 \a 2 n 2 a\; a 2 a'[ ) (15) 

cr i°2 / 

At this point T may be diagonalised to determine A max and hence the M = 2 approximation ( |l2] ) to 
(§■ 

In order to proceed to larger lattices, we must augment (expand) the system and environment 
blocks. We let n < — ► (p, v), p — 1, . . . , to, v =|, \ denote a state for an augmented system block, 
consisting of the initial system block and one site added to the right. There are now m' = 2to = 4 
states: n = 1 < — > (1, f),. . ., n = 4 < — > (2, |). 

The transfer matrix TJ for the augmented system block is defined (see Fig. 4) by 

i s (<7 n \i |er n /j,; a fi ) — / ■i s {o~ P v \o~pv\ a v )T\\y v \fi fi ) (lb) 

v" 

An augmented environment block is defined in a similar way, this time adding a site to the left viz 
Tl{o-'rin'\o-niJL\o-"n") =^T c {a'p'v'\apv:a"v")T2{a"a\v"v) (17) 

Now, the superblock and its associated transfer matrix ([l5]) can once again be formed using the 
augmented blocks with m! 1 — > to, TJ i — > % and 1 — > %. Moreover, the process can be iterated, 
each time augmenting the system and environment blocks to the right and the left respectively. 

However, in order to prevent the superblock basis from becoming too large, we must truncate it 
by capping the number of system (and environment) block states, m, which in principle doubles every 
iteration. To do so, we form reduced density matrices for the augmented blocks by performing an 
appropriate partial trace on the projection operator |^ ma x) (VWxl formed from the eigenstate IV'max) 
of T corresponding to the eigenvalue A max . 



2 



That is, for the augmented system block, the density matrix p s is defined by 
p s (n'\n) = (vi n W2 n 2\ i> max 

){il> 

max 

| a\n\o 2 n 2 ) (18) 

CI U2 

where n < — > (m, 02) and n' < — ► a' 2 ). For the augmented environment block we have 

p c (n'\n) = ^2 {<?ini<7 2 n' 2 \ ip 

max max 

I ain 1 a 2 n 2 ) (19) 

(Tini 

with n < — ► (n 2 ,<J 2 ) and n' < — > (n^cr^). Eigenvalues and eigenvectors of the density matrices are 
then found viz 

{k)W |n) (8) : n= l,...,m'} and K^,|n) (o) : n = l,...,m'} 

for and p c respectively where 1 > > uj^ > ■ • ■ > kfi for i = s or e. 

In proceeding to the next iteration then, we represent the augmented block transfer matrices TJ 
and in terms of the density matrix eigenvectors and truncate the block Hilbert spaces so as to 
retain only the m most important states viz 

m 

(n'\ n"f ) Ti\a , n"n'\ t Tn'"n;a , 'n") (i) (n"'\ n) ^%{a'n' p'\anp;a" p") (20) 

n" ,n"' — 1 

for i = s and e, n, n' = 1, . . . , m and where m is the chosen cap on the number of states per block. 

Now, T is a large, sparse, non-symmetric matrix, for which efficient algorithms for the determina- 
tion of A max and |?/> ma x} are becoming available. So far we have simply used Schur decomposition and 
the inverse power method to perform this task. We have thus only been able to work with m = 16. 
Results for ip, u and Cy for 7 = 1 (pure XY) and 7 = 2 (dimerised XY) are given in Figs 5-7. 

Finite-M results ip^ M ^ are obtained by fixing the inverse temperature (3 = /3q in the expressions (Q) 
and increasing the superblock size 2M, beginning with M — 2, so at each iteration the temperature 
is identified as T = 1/M/3q. For the /3q values we have considered (0.05-0.2) there is a small error due 
to the finiteness of the Trotter decomposition but this is negligible compared with the error due to 
Hilbert space truncation, a result which is easily verified by considering the exact solution for <ip( M ^ 
[0. A typical calculation, using NAG library routines on a 333 MHz DEC Alpha machine, takes 6 
hours to generate a superblock size of 2M = 300. 

We note that convergence at low temperatures is far better for the gapped system (7 = 2). 
However, even in the gapless case, the peak position and peak height in the specific heat are afflicted 
by errors of only ca. 3%, and results may be improved markedly by using larger values of m and the 
finite lattice method jD. Also, results for u and cy are obtained by numerically differentiating spline 
fits of the computed ip^ M ' values, and so improvements may be obtained by combining results from a 
number of different (3q values. 

Note that, at any given iteration, the DMRG usually works with a chain which is spatially finite. 
Here the chain is infinite and the finiteness is in the level of the Trotter approximation. Another 
difference is that the DMRG usually produces its best results for the ground state energy and less 
accurate results for higher excitations. A different situation occurs here — the lower the temperature, 
the less accurate the result. Even so, as can be seen in Fig. 5, the exact ground state energy limx^o ip, 
is recovered with reasonable accuracy, especially in the gapped case 7 = 2. 
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Finally, we note that, as T is non-symmetric, the use of a truncated basis will not in general lead 
to a variational lower bound on A max . However, in practice it can be seen, for instance, that the 
DMRG result for the internal energy u is always bounded above the exact value. 

To conclude, we have applied the DMRG to a quantum spin chain at non-zero temperature, making 
use of Nishimo's adaptation of the method to 2D classical systems. Reasonable results are obtained 
for the specific heat down to low temperatures in calculations involving an extremely small basis 
set, agreement with the exact solution being markedly better in the case where the system has a 
substantial gap. The approach may prove useful in determining thermodynamic properties of models 
of experimentally realisable systems such as coupled chains and models with anisotropy, dimerisation 
and frustration. 

R. J. B. gratefully acknowledges the support of SERC grant no. GR/J26748. Computations were 
performed on the DEC 8400 facility at CLRC. 
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Figure Captions 

1. Pictorial representation of matrices a) Ti^fi 1 . . . fj, 2 ^ 1 ^ 1 . . . a 2M ) and b) ^(fi 1 . . . p? M | cr 1 . . . a 2M ). 
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2. Pictorial representation of the initial system block transfer matrix T s (a'n' /j,'\cr n fi; a"fi"). n and 
n' = 1, . . . , to = 2 are initial and final states of the system block which consists of a single site, a, 
a" and a' = | or j and /x, fi" and // = f or J, are initial, intermediate and final states for the 
adjeacent sites to the left and right of the block respectively. The intermediate state n" of the system 
block is summed over to produce the matrix product. 

3. Pictorial representation of the superblock, which consists of system and environment blocks con- 
nected by two added sites to form a periodic chain. ri\ and and n-i and n' 2 are initial and final 
state indices for the system and environment blocks respectively. <j\ and a[ and a 2 arL d a' 2 are initial 
and final states for the two added sites. The intermediate states a" and a 2 are summed over to form 
the matrix product. 

4. Pictorial representation of the augmented (a) system and (b) environment blocks, which consist of 
the old blocks augmented by adding sites to the right and the left respectively, n and n' = 1, . . . , to' 
are the initial and final states of the new (augmented) block which consists of the old block (p and 
p') augmented with an added site (y and v'). a, a" and a' and /1, and /i' are state indices for the 
adjacent sites. The intermediate state v" of the added site is summed over to produce the matrix 
product. 

5. Exact and DMRG results for the Helmholtz free energy tp as a function of temperature T for two 
cases: 7=1 and 7 = 2. 

6. Exact and DMRG results for the internal energy u as a function of temperature T for two cases: 
7=1 and 7 = 2. 

7. Exact and DMRG results for the specific heat cy as a function of temperature T for two cases: 
7=1 and 7 = 2. 
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